First order phase transitions in classical lattice gas spin models 



H. Chamati 

Institute of Solid State Physics, Bulgarian Academy of Sciences, 
72 Tzarigradsko Chaussee, 1784 Sofia, Bulgaria. 



o 
o 

(N 



in 



o 

S 



I 

o 



> 



(N 

in 
o 
l> 
o 



X 



S. Romano 

?7njfa di Ricerca CNISM e Dipartimento di Fisica "A. Volta", 
Universitd di Pavia via A. Bassi 6, 1-27100 Pavia, ITALY 
(Dated:) 

The present paper considers some classical ferromagnetic lattice-gas models, consisting of particles that 
carry n-component spins (n = 2, 3) and associated with a _D-dimensional lattice (D = 2, 3); each site can 
host one particle at most, thus implicitly allowing for hard-core repulsion; the pair interaction, restricted to 
nearest neighbors, is ferromagnetic, and site occupation is also controlled by the chemical potential ji. The 
models had previously been investigated by Mean Field and Two-Site Cluster treatments (when D — 3), as 
well as Grand-Canonical Monte Carlo simulation in the case /i = 0, for both D — 2 and D = 3; the obtained 
results showed the same kind of critical behaviour as the one known for their saturated lattice counterparts, 
corresponding to one particle per site. Here we addressed by Grand-Canonical Monte Carlo simulation the case 
where the chemical potential is negative and sufficiently large in magnitude; the value fi = —D/2 was chosen 
for each of the four previously investigated counterparts, together with /i — — 3D/4 in an additional instance. 
We mostly found evidence of first order transitions, both for D = 2 and D = 3, and quantitatively characterized 
their behaviour. Comparisons are also made with recent experimental results. 

PACS numbers: 75.10.Hk, 05.50.+q, 64.60.-i 



I. INTRODUCTION 

There exist a few statistical mechanical models involving 
classical continuous spins, by now extensively investigated, 
and which play a central role in a variety of real physical sit- 
uations; they have been especially studied in their saturated 
lattice (SL) version, where each lattice site hosts one spin. 
The interest in their lattice gas (LG) extensions was recently 
revived, and they were addressed by means of analytical the- 
ories and simulation (see for example HH^^IHHSSSIMl 
and references therein). LG spin models are obtained from 
their SL counterparts by allowing for fluctuations of occu- 
pation numbers, also controlled by the chemical potential fi. 
These models have often been used in connection with alloys 
and absorption; the methodology somehow allows for pres- 
sure and density effects. 

As for symbols and definitions, classical SL spin mod- 
els involve component unit vectors u^, associated with a 
D— dimensional (bipartite) lattice Z^; let denote dimen- 
sionless coordinates of the lattice sites, and let Uk,a denote 
cartesian spin components with respect to an orthonormal ba- 
sis Ba, whose unit vectors can be taken as defined by the 
lattice axes. The orientations of the magnetic moments of 
the particles are parameterized by usual polar angles {(f>j} 
(n = 2) or spherical ones {{(pj, 0j)} {n ~ 3). The inter- 
action potential, restricted to nearest neighbors, is assumed to 
be ferromagnetic and, in general, anisotropic in spin space, 
i.e. 



jk 



<^Qjk, Qjk 



Uk., 



e > 0, a > 0, 6 > 0, a + b> 0, max(a, b) = 1. (1) 
Notice also that the condition max(a, &) = 1 in the above 



equation can always be satisfied by a suitable rescaling of e; 
here and in the following the quantity e will be used to set tem- 
perature and energy scales; thus T ~ kst/e, where t denotes 
the absolute temperature and fcs is the Boltzmann constant; 
the corresponding (scaled) Hamiltonian is given by: 



{j<h} 



(2) 



The case n — 1 corresponds to the Ising model; isotropic 
0(n)-syinmetric models (n > 1) correspond to a = b, 
Qjk ~ —Uj ■ Uk, and are referred to as planar rotators (PR, 
n = 2) or classical Heisenberg model (He, n = 3); the ex- 
tremely anisotropic and 0(2)-symmetric XY model is defined 
by ?i = 3, a = 0. For these models the simplification resulting 
from the neglect of translational degrees of freedom makes it 
possible to obtain rigorous mathematical results lol flolfTlll en- 
tailing existence or absence of a phase transition, and, on the 
other hand, to study it by a whole range of techniques, such 
as Mean Field (MF) and Cluster Mean Field treatments, high- 
temperature series expansion of the partition function, Renor- 
malization Group (for a recent review see lfl2|] ). computer 
simulation (usually via Monte Carlo (MC) methods ifisll ). 

LG extensions of the continuous-spin potential model con- 
sidered here are defined by Hamiltonians 



{j<k} 



N 



E 



Vk, 



(3) 



where Vk = 0,1 denotes occupation numbers; notice that 
A < reinforces the orientation-dependent term, whereas 
A > opposes it, and that a finite value of A only becomes 
immaterial in the SL limit /i +oo. It is worth mentioning 
that in such systems the fluctuating occupation numbers give 
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rise to additional fluid-like observables in comparison to the 
usual SL situation. 

Rigorous results entailing existence or absence of an or- 
dering transition are also known for LG models with contin- 
uous spms H III H [B H. For some models defined 
hy D ~ 3, interactions isotropic in spin space, and sup- 
porting a ferromagnetic phase transition in their SL version, 
it has been proven that there exists a /ig, such that, for all 
/i > /io- the system supports a ferromagnetic transition, with 
a /^t— dependent transition temperature. Notice that /io < 
when A < lfT5l[l6l[T7ll . whereas a positive /^to may be needed 
when A > 0. More recently ifisll . the existence of a first-order 
transition, involving discontinuities in both density and mag- 
netization, has been proven for the isotropic case (and D = 3), 
in a suitable regime of low temperature and negative fi. 

For D = 2, the SL-PR model produces at low- 
temperature the extensively studied Berezinskii-Kosterlitz- 
Thouless (BKT) transition lfl9i I20II : the existence of such a 
transition for the LG counterpart has been proven rigorously 
as well ll2lll . More recently, it was rigorously proven Il22ll that, 
for /i negative and sufficiently large in magnitude, the transi- 
tion becomes first-order 

Notice also that the above mathematical theorems do not 
yield useful numerical estimates of the ji value where the 
change of transition sets in; some answer to this question can 
be looked for by analytical approximations such as MF or 
Two-Side-Cluster (TSC) treatments mUt], or by simulation 

aai. 

The Hamiltonian (Eq. (|3]l) can be interpreted as describ- 
ing a two-component system consisting of interconverting 
"real" (I'k = 1) and "ghost", "virtual" or ideal-gas parti- 
cles (h'k = 0); both kinds of particles have the same ki- 
netic energy, fj, denotes the excess chemical potential of "real" 
particles over "ideal" ones, and the total number of parti- 
cles equals the number of available lattice sites (semi-Grand- 
Canonical interpretation). The semi-Grand-Canonical inter- 
pretation was also used in early studies of the phase dia- 
gram of the two-dimensional planar rotator, carried out by 
the Migdal-Kadanoff RG techniques, and aiming at two- 
dimensional mixtures of '^He and ''He ll23i l24ll . where non- 
magnetic impurities correspond to '^He. 

In the three-dimensional case, the topology of the phase di- 
agram of the model (O had been investigated by MF and TSC 
approximations for the Ising 101 as well as PR cases ^ in the 
presence of a magnetic field, and for He at zero magnetic field 
lit]. These investigations were later extended |f^ to extremely 
anisotropic (Ising-like) two-dimensional LG models defined 
by a = 1, 5 = in Eq. ([T), and in the absence of a mag- 
netic field as well. The studied models were found to exhibit 
a tricritical behaviour i.e. the ordering transition turned out 
to be of first order for fi below an appropriate threshold, and 
of second order above it. When the transition is of first or- 
der, the orientationally ordered phase is also denser than the 
disordered one. For the three-dimensional PR these finding 
were confirmed, recently, by simulation in connection with the 
phase diagram of He It has been found that, despite the 
simplicity of LG spin models, their predictions broadly agree 
with the ones obtained by means of more elaborate magnetic 



fluid models (see e.g. ||25|] and references therein). 

On the other hand, thermodynamic and structural properties 
had been investigated by means of Grand-Canonical Monte 
Carlo simulations as well ^ tsj], for particular values of the 
chemical potential equal or close to zero. It had been found 
that there is a second order ferromagnetic phase transition 
manifested by a significant growth of magnetic and density 
fluctuations. The transition temperatures were found to be 
about 20% lower than that of the corresponding SL values and 
the critical behaviour of the investigated models to be consis- 
tent with that of their SL counterparts. Furthermore it had 
been found that MF yields a qualitatively correct picture, and 
the quantitative agreement with simulation could be improved 
by TSC, which has the advantage of predicting two-site cor- 
relations. 

Notice also that the above Hamiltonian (Eq. (|3]l) describes 
a situation of annealed dilution; on the other hand, two- 
dimensional models in the presence of quenched dilution, and 
hence the effect of disorder on the BKT transition, have been 
investigated using the PR model Jgl |27i |28,. ,29, ,30l Isill and 
very recently its X Y counterpart Bill ; it was found that a suf- 
ficiently weak disorder does not destroy the transition, which 
survives up to a concentration of vacancies close to the per- 
colation threshold. Let us also remark that two-component 
spins are involved in the PR case, whereas XY involves three- 
component spins but only two of their components are in- 
volved in the interaction: in this sense the two models en- 
tail different anchorings with respect to the horizontal plane 
in spin space. Two-dimensional annealed lattice models were 
investigated ||8{| as well, and the obtained results for /i = or 
a moderately negative /i were found to support those obtained 
for quenched models. For a large negative /i, renormalization 
group treatments had suggested ll23l[24ll that the transition be- 
tween the BKT and the paramagnetic phase is of first order 

In this paper, we present an extensive Monte Carlo study 
of some LG ferromagnetic models, where ji is negative and 
comparatively large in magnitude (notice that ji < —D would 
produce an empty ground-state), in order to gain insights into 
their critical behaviour and to check the impact of the chemi- 
cal potential on their physical properties. On the other hand, 
for D = 3, we will also test the MF or TSC approximations 
used to obtain the phase diagrams of Refs. yl tSfl. In keeping 
with our previous studies, the models are further simplified 
by choosing A = 0, i.e. no pure positional interactions. As 
for the values of the chemical potential, we chose fi = —D/2 
in the four cases corresponding to our previous investigations 
with /i = 0, and carried out additional simulations for D — 2, 
PR and ^ = — 3Z)/4, as explained below. 

The rest of the paper is organized as follows; the simulation 
procedure is briefly explained in Section [III section |III] ana- 
lyzes the simulation results. Finally, the effects caused by the 
chemical potential on the nature of the transition are discussed 
in SectionHVl which summarizes our results, and where some 
comparisons are made with a recent experimental work. 
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n. MONTE CARLO SIMULATIONS 

A detailed treatment of Grand-Canonical simulations can 
be found in or via Refs. flll5[ [32ll : the method outlined here 
has already been used in our previous studies of other LG 
models IH Is], |^. Simulations were carried out on periodi- 
cally repeated samples, consisting of V = sites, where 
L = 40, 80, 120, 160 for D = 2, and L = 10, 20, 30 for 
D ~ 3, i.e. in keeping with the named previous studies of 
ours; calculations were carried out in cascade, in order of in- 
creasing reduced temperature T. 

The two basic MC steps used here were Canonical and 
semi-Grand-Canonical attempts; in addition two other fea- 
tures were implemented II33L l34jl : (i) when a lattice site was 
visited, Canonical or semi-Grand-Canonical steps were ran- 
domly chosen with probabilities Vcan and Vgc, respectively; 
we used T'can/^GC = ?i — 1, since spin orientation is defined 
by (n — 1) angles, versus one occupation number and (ii) sub- 
lattice sweeps (checkerboard decomposition) lf33[ [3; thus 
each sweep (or cycle) consisted of 2V attempts, first V at- 
tempts where the lattice sites was chosen randomly, then V/2 
attempts on lattice sites of odd parity, and finally V/ 2 attempts 
on lattice sites of even parity. Equilibration runs took between 
25 000 and 200 000 cycles, and production runs took between 
250 000 and 1 000 000; macrostep averages for evaluating sta- 
tistical errors were taken over 1 000 cycles. Different random- 
number generators were used, as discussed in Ref. lf34ll . 

Computed thermodynamic observables included 
mean Hamiltonian per site, H = (1/F)(A), den- 
sity p = i^/V) (N), as well as their derivatives 
with respect to temperature or chemical potential, 

C^v/kB = {l/V){d{A)/dT)^y, pT = [dp/dT)^y, 

Pfi = {dp/ diJi)Ty, defined by appropriate fluctuation 
formulae 113211 . 

We also calculated mean magnetic moment per site and sus- 
ceptibility, defined by 

7\/ ^ 1 (yr^'^ , (4) 

where for PR or He the vector F is defined by 

V 

F = ^z/feUfc, (5) 

k=l 

whereas only the in-plane components of the vector spins (i.e. 
only the Cartesian components explicitly involved in the inter- 
action potential) are accounted for in the XY case. 

The behaviour of the susceptibility was investigated by con- 
sidering the two quantities: 

Xi = ^((F.F)-(|F|)^) (6) 

and 

X2 = 4 (F ■ F) ; (7) 



simulation estimates of the susceptibiUty ifTsi [35i [36|] are de- 
fined by 

_ J Xij in the ordered region 
^ 1x2, in the disordered region ' 

notice also that, for a finite sample, X2 < PV, and that x = X2 
in two dimensional cases. 

A sample of V sites contains DV distinct nearest- 
neighbouring pairs of lattice sites; we worked out pair oc- 
cupation probabilities, i.e. the mean fractions Rjk of pairs 
being both empty (Ree = ((1 — ~ Vk))), both occupied 
{Roo = {vjVk)), or consisting of an empty and an occupied 
site {Reo = ((1 — Vj)vk + (1 — Vk)vj)). It should be noted 

that i?ee + Roo + Reo = 1- 

Short- and long-range positional correlations were com- 
pared by means of the excess quantities 

<o - In [j^^ , Ko = Roo - p\ (9) 

collectively denoted by R*^ (notice that these two definitions 
entail comparable numerical values). 

Quantities such as p, px, and the above pair correla- 
tions Rjk or R*g can be defined as "fluid-like", in the sense 
that they all go over the trivial constants in the SL limit. Let 
us also remark that some of the above definitions (e.g. 
C^.y and pt) involve the total potential energy both in the 
stochastic variable and in the probability measure ("explicit" 
dependence), whereas some other definitions, e.g. p^ or the 
quantities Rjk, involve the total potential energy only in the 
probability measure ("implicit" dependence). 

m. SIMULATION RESULTS 

A. D = 2,PR,p = -1 

Simulation results, obtained in the named cases for a num- 
ber of observables, such as the mean energy per site and den- 
sity, showed that these quantities evolve with the tempera- 
ture in a smooth way, and were found to be independent of 
sample sizes. In comparison with Ref. fst], their temperature 
derivatives C^y and pT (Fig. [T]i showed recognizably more 
pronounced peaks about the same temperature T « 0.51, 
around which the sample size dependence of results became 
slightly more pronounced. Comparison with our previous re- 
sults shows that the location of the maximum of C^v is shifted 
towards lower temperatures as p < grows in magnitude. 

Plots of lnx2 versus T, reported in Fig. |2] show results 
independent of sample size for T > 0.52, and then their pro- 
nounced increase with sample size for T < 0.51, suggesting 
its divergence with L. In general this case qualitatively re- 
produces our previous simulation results 180, but with more 
pronounced derivatives and peaks at a lower temperature. In 
order to estimate the critical temperature we applied the finite 
size scaling analysis, along the lines discussed in reference 
lit], and here again we found a BKT transition occurring at 
Tbkt = 0.502 ± 0.002, corresponding to a particle density 
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about 0.832 ± 0.003. Comparison between our previous re- 
sults and the present ones shows that both transition temper- 
ature and "critical" particle concentration are monotonically 
increasing with the chemical potential (see Table|I]i. 
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FIG. I: Simulation estimates for the specific heat per site C^v ver- 
sus temperature, obtained with different sample sizes, for the two- 
dimensional LG-PR and = — 1. Simulation results for pr and 
obtained with the largest examined sample size are shown on the top. 
Statistical errors range between 1 % and 5 %. 
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FIG. 2: Simulation estimates for the logarithm of the magnetic sus- 
ceptibility X2 versus temperature, obtained with different sample 
sizes, for the two-dimensional LG-PR and /i = — 1. Unless oth- 
erwise stated, here and in the following figures, statistical errors fall 
within symbol sizes. 



rapid variations across the transition region, in accordance 
with the sharper maximum of p^. Their behaviours suggest 
inflection points roughly corresponding to the maximum of 

Short- and long-range positional correlations have been 
compared via the excess quantities i?*^, whose simulation re- 
sults for the largest sample size are shown in Fig. [H show- 
ing sharper maxima than their counterparts corresponding to 
pt = 0. Notice also that the position of the maximum for i?"^ 
again corresponds to the location of the peak of C^v. The 
quantities i?*^ are rather small, and this could be traced back 
to the absence of pure positional interactions. 



For the SL-PR model the maximum of the specific heat is 
located at about 15% ll37l]a bove the BKT transition; for the 
LG-PR model and p = | s]] we had found a broad peak about 
5% above the BKT transition, and here we find a sharper one 
about 2% above the transition temperature. 

For p = — 1, fluidlike quantities show qualitatively similar 
behaviours as their counterparts obtained for pt = 0. Results 
for Pt and p^, obtained with the largest sample sizes are re- 
ported in Fig. [T] they were found to behave in a similar fashion 
to the specific heat and to exhibit sharper peaks taking place 
at the same temperature as that of C^y. Recall that p^ has a 
broad maximum for p = 0. In other words, here the ferro- 
magnetic orientational fluctuations taking place in the transi- 
tion range do produce stronger fluctuations of site occupation 
variables, and this tends to reduce the difference between "im- 
plicit" and "explicit" dependencies on the potential energy as 
mentioned in Ref . . 

Pair occupation probabilities Rjk were found to be insen- 
sitive to sample sizes; results for our largest sample size are 
shown in Fig. |3] These quantities are monotonic functions of 
temperature as their counterparts for p = 0, but with more 
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FIG. 3: Simulation estimates for the pair occupation probabilities 
Rjk versus temperature, for the two-dimensional LG-PR with lin- 
ear sample size L = 160. The results refer to p = —1. 
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TABLE I: Transition temperatures and "critical" particle density pc of PR and XY models for some selected values of the chemical potential 
fj,. Depending on fj, we have either a BKT or a first order one (I); here pc denotes the density at the BKT transition temperature. 



Model 


u 


Transition 


e 


0^ 
r c 




OO 


13 jV 1 




1 

1 . 




0.1 


BKT 


0.75 ±0.01[a] 


0.938 ± 0.002 




0.0 


BKT 


0.733 ± 0.003[8J 


0.924 ± 0.003[8] 




-0.2 


BKT 


0.71 ±0.01[a] 


0.900 ± 0.002 




-1.0 


BKT 


0.502 ± 0.002 


0.832 ± 0.003 




-1.5 


I 


0.279 ± 0.001 




XY(n=3) 


OO 


BKT 


0.700 ± 0.005|31] 


1. 




0.0 


BKT 


0.574 ± 0.003[8J 


0.918 ± 0.004[8] 




-1.0 


I 


0.333 ±0.001 
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FIG. 4: Simulation estimates for the excess quantities Roo for the 
two-dimensional LG-PR. Simulation results were obtained with 
L = 160 and = -1. 



B. D — 2 and first-order transitions 

Additional simulations carried out for D = 2, PR, i^i = 
—3D /A, showed a recognizably different scenario. Here, for 
all investigated sample sizes, we found pronounced jumps 
of various observables, such as H, p (FIG. H) and even M 
(which kept decreasing with increasing sample size), taking 
place over a narrow temperature range, AT = 0.0005. No- 
tice that X2 remains independent of sample sizes in the high- 
temperature regime, and then develops a pronounced increase 
with sample size. From a comparison of the behaviours of X2 
for /i = — 1 (Fig. [2]) and ji = —1.5 (Fig. |7]) one can ob- 
serve the change of the critical behaviour at the two values 
of n. For /i — —1.5 the thermodynamic observables show a 
discontinuous behaviour characteristic of a first-order transi- 
tion, now to a low-temperature BKT phase. The behaviours of 
Cf_iVi Pt, Pfj. 3re shown in Fig. |6] and also exhibit pronounced 
differences from their counterparts in the previous case (see 
also below). 

This result confirms previous RG predictions lUEHSl; on 
the other hand, recent simulation studies addressing quenched 
dilution have found that the transition temperature vanishes 



below the percolation threshold S HH HI [H . 

Notice that usage of the Grand-Canonical ensemble allows 
quite wide changes of density with temperature; in the inves- 
tigated cases we used ji > —D, and found that p « 1 in the 
low-T phase, where pT < 0; such changes are obviously ex- 
cluded from the start in the treatment of a quenched-dilution 
model. On the other hand, values p, < —D produce an essen- 
tially empty ground-state; in this regime one can expect that 
p to increase with T, only becoming appreciable above some 
threshold, and that the BKT phase disappears. 
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FIG. 5: Simulation results for the density p and the mean energy 
per site —H obtained for the two-dimensional LG-PR. The value 
/i = —1.5 was used for the present simulations. 

Here and in the following subsection, transitional proper- 
ties such as AH, Ap, (as well as A/, in the next subsection), 
were estimated by analyzing simulation results for the largest 
sample size as discussed in Refs. llssilsoll . The relevant results 
are reported in Table HIl 

Let us now turn to the discussion of the nature of the low- 
temperature phase. Here the magnetization was found to ex- 
hibit a power-law decay with increasing sample size. A fit to 
the expression 

lni\/ = -6ilnL + 6o, h > 0: (10) 

showed that the ratio hi{T)/T is a constant. This shows that 
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FIG. 6: Simulation estimates for tiie specific lieat per site C^v ver- 
sus temperature, obtained with different sample sizes for the two- 
dimensional LG-PR and /i — —1.5. Simulation results for pr and 
Pfi obtained with the largest examined sample size are shown on the 
top. Statistical errors range between 1 % and 5 %. 
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FIG. 7: Simulation estimates for the logarithm of the magnetic sus- 
ceptibility X2 against temperature, obtained with different sample 
sizes for the two-dimensional LG-PR and p = —1.5. 



the magnetization goes to zero in the thermodynamic limit 
(L oo), as predicted by the Mermin Wagner theorem for 
2D systems, where no long range order should survive. Note 
that this behaviour is consistent with the spin wave theory 
developed for the two-dimensional saturated planar rotator 



TABLE II: A summary of simulation estimates for properties at first- 
order transition for the two-dimensional models obtained using sim- 
ulations. 



Model 




e 


AH 


Ap 


PR 


-1.5 


0.279 ± 0.001 


0.3562 ± 0.0005 


0.9917 ±0.0001 


XY 


-1.0 


0.332 ± 0.001 


0.664 ±0.002 


0.910 ±0.001 



model iHIMl. 

Results for In X2 against temperature (Fig. |7]) were found to 
be independent of sample size when T > 0.281, and showed 
a recognizable increase with it (a linear dependence of In X2 
on InL) when T < 0.278. Thus in the low-temperature re- 
gion the susceptibility exhibits a power law divergence with 
the linear sample size, showing a BKT phase ll37ll40ll . 

As for simulation results obtained for the XY LG model 
with fi = —1, it was found that the thermodynamic quantities 
have qualitatively similar behaviours as those obtained for the 
above LG-PR with /i = —1.5. The phase transition was found 
to be first order taking place at T = 0.332 ± 0.001; estimates 
of transition temperatures reported in Table |T] show that they 
increase as a function of the chemical potential. Transitional 
properties of the mean energy, the density and the magnetiza- 
tion are presented in Table HIl Here again we have found that 
the transition takes place from a paramagnetic to a BKT-like 
phase. 

Let us emphasize that, as remarked above, PR and XY mod- 
els entail different anchorings with respect to the horizontal 
plane in spin space; this difference correlates with the pro- 
nounced qualitative different in transition behaviour observed 
when /i = — 1 . 

When both PR and XY lattice gas models exhibited a first 
order phase transition, their fluidlike quantities were found to 
behave in a qualitatively similar way. The following discus- 
sion will concentrate on these properties for the XY model. 

Fig. |6] shows that pr is negative and decreases with in- 
creasing T in the low-temperature region (where it is essen- 
tially driven by orientational correlations), and then it be- 
comes weakly positive and increasing with T in the high- 
temperature phase; thus, here and in the following subsection, 
p decreases with T in the low-temperature phase, and then 
increases with T in the high-temperature region. On the other 
hand, here is an increasing function of T, exhibiting a jump 
across the transition. 

Simulation results for the pair occupation probabilities, re- 
ported in Fig. [8] and the excess quantities i?*^ shown in Fig. 
|9] reveal that these quantities are discontinuous at the first or- 
der transition region. On the other hand they show the effects 
caused by the ferromagnetic interaction on the density in the 
system. The quantity i?"^ remains negligible due to the ab- 
sence of purely positional interaction. The behaviour of these 
quantities follow in general the trends of the mean Hamilto- 
nian and the density. To summarize we found that the system 
exhibits a first order phase transition form a dense BKT phase 
to a paramagnetic one; in the temperature-density phase dia- 
gram, both phases are expected to coexist over some range of 
densities and temperatures. 
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FIG. 8: Simulation estimates for thie thuee pair occupation probabili- 
ties RjK- for the two-dimensional LG-XY, for a sample with linear 
size L = 160. The value = — 1 was used in this simulation. 
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TABLE III: Transition temperatures B and "critical" particle density 
Pc of PR and He models for some selected values of the chemical 
potential /j.. Depending on fi, there is either a second order transition 
(II) or a first order one (I); pc denotes the density at the second order 
transition temperature. 



Model 




Transition 


e 


Pc 


PR(n=2) 


oo 


II 


2.201 ±0.003 


1. 




0.1 


II 


1.423 ±0.003 


0.6900 ± 0.004 




-1.5 


I 


0.794 ±0.001 




He(n=3) 


oo 


II 


1.443 ±0.001[36] 


1. 




0.0 


II 


0.998 ± 0.001[5] 


0.743 ± 0.002 




-1.5 


I 


0.557 ±0.001 





describes qualitatively well the transitional properties of the 
named models and that TSC improves upon it. In Tablelllllwe 
report results for the transition temperature obtained, via sim- 
ulations, for some values of fi for these models so far Here 
also one can read that the transition temperature decreases 
with decreasing fi. 

Simulation results for both models exhibited a recognizable 
qualitative similarity, so that only plots of PR are presented 
here. Behaviours of observables such as mean energy, density 
p and magnetisation M (shown in Fig. [TOl l were found to 
be either size independent or to depend slightly on sample 
sizes in the transition region. Furthermore, for all examined 
sample sizes, we found abrupt jumps of these observables, 
taking place over a narrow temperature range, AT = 0.0005. 
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FIG. 9: Simulation estimates for the quantities Roo, obtained with 
L = 160 for the two-dimensional LG-XY and p. — —1. 
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C. D — 3 and first-order transitions 



Simulation results presented in this subsection for the 
three-dimensional PR and He models show the effects caused 
by large negative pt on their transitional behaviour, and, on 
the other hand, can be used to check the predictions of the 
molecular-field like treatments used to construct the phase di- 
agrams reported in our previous papers iHHt] ; we refer to them 
for further details, and present here only the final numerical 
results for the specific cases of interest. 

It is well known that these approximate treatments do not 
describe fluctuations adequately, so that their predictions have 
to be taken with caution. For example, MF predicted a first 
order phase transition at ^ = 0, while TSC and MC gave 
evidence of a second order phase transition for He iU. For 
both three-dimensional models, simulations performed for a 
selected value of the chemical potential, revealed that MF 



0.0 I • 1 • ^^^^"^WfWYwW M « 

0.77 0.78 0.79 0.80 0.81 

T 

FIG. 10: Simulation results for the magnetisation M against the 
temperature for the three-dimensional LG-PR, obtained with differ- 
ent sample sizes and ^ = — 1.5. 

In Table |IV] we present the transitional properties such as 
jumps in mean energy per site and density, respectively, as 
well as magnetisation in the ordered phase, at the first order 
phase transition undergone by the three-dimensional PR and 
He; these results were obtained via MC, MF and TSC. Com- 
parison shows that TSC produces a better estimate than MF 
for the transition temperature; on the other hand, MF better 
predicts the jumps of thermodynamic quantities at the transi- 
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tion. In general, according to the results gathered in Table Hill 
one can see that the phase diagram predicted by the approxi- 
mate molecular field theories is at least qualitatively correct. 
This fact is confirmed by the recent simulation results for the 
phase diagram of the diluted PR reported in Ref . 0] . 

The susceptibility, actually xi, reported in Fig. (TT] showed 
a peak at a temperature about 0.792, a strong sample size 
dependence below this temperature and no sensitivity to the 
sample sizes above it. The behaviours of the three derivatives 
C^v, PTi Pfj. (not reported) were found to be qualitatively 
similar to Fig. |6]l. 




0.81 



FIG. 11: Simulation results for the susceptibility xi for the three- 
dimensional LG-PR, obtained with different sample sizes. The as- 
sociated statistical errors, not shown, range up to 10%. The value 
fi = —1.5 was used in this simulation. 

Other fluidlike quantities such as Rjk (Fig [T2] ) and i?*^ 
(not reported here) show how the density behaves when the 
three-dimensional PR lattice gas model exhibits a first order 
transition. These quantities are discontinuous at the transition 
temperature and follow the behaviour obtained for the den- 
sity and the mean energy; once more we witnessed the small- 
ness of the excess quantities Roo due to the absence of purely 
positional interaction. In general we remarked a pronounced 
qualitative similarity between the behaviours of the fluidlike 
quantities in the present case and those discussed in the previ- 
ous subsection for 2D models. 



IV. CONCLUDING REMARKS 

We have studied the critical properties of four LG models 
defined hy fi < and sufficiently large in magnitude, fj, = 
-D/2, plus an additional case (D = 2, PR, ^ = -3L'/4). 
This allowed us to investigate the impact of the chemical po- 
tential on the nature of phase transition of these models and 
thus to gain insights into their phase diagrams. Our simu- 
lations were performed in the absence of pure positional in- 
teraction. A number of thermodynamic quantities including 
some characteristics of fluid systems were estimated. It was 
found that the common feature of most cases is the onset of a 




0.81 



FIG. 12: Simulation results for the three pair occupation probabil- 
ities Rjk obtained for the three-dimensional LG-PR with linear 
sample size L = 30. The value jj, = —1.5 was used in this sim- 
ulation. 



first-order phase transition induced by the ferromagnetic in- 
teraction, and where an an abrupt change in the density of the 
system was observed. 

In two-dimensions we have investigated both PR and XY 
models for n = —D/2. At this value of the chemical poten- 
tial they showed different critical behaviours. PR exhibited 
a BKT phase transition, while XY showed a first order one. 
This might be a consequence of the fact that the two models 
entail different anchorings with respect to the horizontal plane 
in spin space. PR was further studied for fi = — 3Z5/4, where 
evidence of a first order transition was found. The change of 
the nature of the phase transition from BKT to a discontinu- 
ous one agrees with previous RG predictions ll23l l24ll and rig- 
orous mathematical results [2Z]; on the other hand, in recent 
simulation studies of quenched dilution it was found that the 
transition temperature vanishes below the percolation thresh- 
old S HI, US IH- 

Notice that usage of the Grand-Canonical ensemble allows 
quite wide changes of density with temperature. Such changes 
are obviously excluded from the start in the treatment of a 
quenched-dilution model. Thus, there are significant differ- 
ences between both methods, yet the two resulting pictures 
are somehow compatible. 

Phase transition and critical dynamics in site-diluted arrays 
of Josephson junctions were recently studied experimentally 
in Ref. ll4lll : according to the Authors' results, the BKT tran- 
sition is altered by the introduction of percolative disorder far 
below the percolation threshold. Furthermore, the Authors of 
Ref. ll4lll found evidences of a non-BKT-type superconduct- 
ing transition for strongly disordered samples, taking place at 
finite temperature. Our results suggest that the transition in 
the named region becomes of first order 

For the three dimensional models investigated here, i.e. PR 
and He, we found a first order phase transition form a fer- 
romagnetic dense phase to a diluted paramagnetic one. The 
results obtained via simulation for /i = —D/2 were found 
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TABLE IV: Estimates for some properties at first-order transition for tlie tliree-dimensional PR and He obtained by different approaches. The 
results are obtained with ^ — —1.5. 

Model Method 6 AH Ap M 

~PR MC 0.794 ±0.001 0.910 ± 0.004 0.684 ± 0.002 0.772 ± 0.002 

MP 0.741 1.138 0.849 0.897 

TSC 0.760 1.518 0.756 0.903 

He MC 0.557 ±0.001 0.882 ± 0.003 0.877 ± 0.001 0.804 ± 0.001 

MP 0.462 0.944 0.958 0.9126 

TSC 0.482 0.959 0.786 0.888 



to confirm those obtained by molecular field aaproximations 
used to construct the phase diagram of Refs. yllSl], showing 
that the phase diagrams obtained there are qualitatively cor- 
rect. 
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